Spin susceptibility of underdoped cuprates: the case of Ortho-II YBa2Cu306.5 
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Recent inelastic neutron scattering measurements found that the spin susceptibihty of detwinned 
and highly ordered ortho-II YBa2Cu306.5 exhibits, in both the normal and superconducting states, 
one- dimensional incommensurate modulations at low energies which were interpreted as a signature 
of dynamic stripes. We propose an alternative model based on quasiparticle transitions between the 
arcs of a truncated Fermi surface. Such transitions are resonantly enhanced by scattering to the 
triplet spin resonance. We show that the anisotropy in the experimental spin response is consistent 
with this model if the gap at the saddle points is anisotropic. 

PACS numbers: 74.72-h,75.25+z,76.50+g 



Recent experiments^ '^'^ have addressed the magnetic 
spectrum in detwinned and highly ordered ortho-II 
YBa2CusOQ,^ which has an average doping per planar 
Cu X ~ 0.09 and Tc ~ 60 K. The ortho-II phase is charac- 
terized by a periodic alternation of filled and empty Cu- 
O 6-axis chains, doubling the size of the unit cell in the 
a direction^. The oxygen ordering reduces the disorder 
making it an ideal candidate for the study of the magnetic 
properties at underdoping. The magnetic response shows 
a resonance^ at 33 meV at Q = (tt, tt) and a rm^-like 
high energy branch dispersing upwards^. At lower en- 
ergies, there is inelastic scattering from one-dimensional 
modulations at incommensurate wavevectors stretching 
down to zero energy but without elastic or quasielastic 
peaks. On entering the superconducting state, the scat- 
tering at energies less than 16 meV is suppressed, but 
not eliminated. In this strongly underdoped regime the 
incommensurability 5 at low energies is very small. Con- 
trary to what has been observed at optimal doping the 
incommensurability below the resonance persists in the 
normal pseudogap phasA The feature associated with 
the one-dimensional modulations, which were interpreted 
as a signature of dynamic stripes^, is a flat top lineshape 
centered at (tt, tt) along the direction perpendicular to 
the chains. An anisotropic response, but no sign of one- 
dimensionality has been reported by Hinkov et on 
untwinned YBa2Cu306+ic. 

Static stripes should show up as elastic magnetic peaks 
and quasi-elastic peaks would be expected^, from slow 
fluctuating stripes, but these are not observed. Also, the 
dramatic drop of the magnetic scattering at low energies, 
at the onset of superconductivity, is not expected if or- 
dered antiferromagnetic regions between hole stripes are 
involved. The stripe-based model also has difficulties to 
explain the isotropy found above the resonancA 

Many authors have used RPA^ (random phase approx- 
imation) schemes to describe the magnetic response near 
optimal doping. This model relates the resonance peak 
and inconmensurate branches with a particle-hole collec- 
tive mode bounded below the superconducting gap. An 
extension to the normal pseudogap phase at underdop- 
ing can be made by assuming an energy gap from an- 



other origin (e.g. d-density wave order^) but even then 
careful tuning of the parameters is required since the res- 
onance decreases while the gap increases in energy with 
decreasing doping. Alternatively one can interpret the 
resonance as a collective mode of a spin liquid which 
goes soft as the antiferromagnetic instability at strong 
underdoping is approachedifl or as an undamped spin- 
wave excitation brought about by strong antiferromag- 
netic correlations^^. 

In this letter we introduce a simple phenomenologi- 
cal model for the pseudogap phase and focus on the 
anisotropic incommensurate low energy response. We 
start from a Fermi surface of quasiparticles which is trun- 
cated near the saddle points leaving only four Fermi arcs 
centered on the diagonals, similar to the form observed in 
ARPES^^ experiments. Although the dispersion we use 
is derived from a staggered flux phase^^-^^, we propose 
it as a phenomenological form for the normal state of a 
doped spin liquid (SL). Arguments that a Fermi surface 
truncation can occur without symmetry breaking have 
been given by Ledermann et a!^^, based on an analysis 
of multi-leg Hubbard ladders and by Honerkamp et 
and also Lauchli et al^^ from a functional renormaliza- 
tion group analysis of a 2-dimensional Hubbard model. 
In addition to the quasiparticle scattering we introduce a 
resonant spin triplet mode at (tt, tt) dispersing upwards, 
in both normal and superconducting phases. We view 
the resonance as the collective mode of the spin liquid 
analogous to that found, for example, in coupled ladder 
systemsiSiifiiSft and in the slave boson approachifl. 

We propose that the low-energy spin response is due 
to particle-hole excitations across the Fermi arcs. In the 
absence of other effects the bare spin susceptibility, due 
to particle-hole excitations, is too small to explain the 
experimental data and has a much richer structure than 
the one observed. In our model such transitions are reso- 
nantly enhanced by scattering into the triplet spin reso- 
nance. The enhancement is strongest close to (7r,7r). We 
analyze the possible sources of anisotropy between direc- 
tions parallel and perpendicular to the chains and flnd 
that results similar to those found experimentally can be 
obtained if the gap at the saddle points is sufficiently 
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anisotropic. 

We start from a two-dimensional dispersion with hop- 
ping to first and second nearest neighbors. Truncation is 
introduced through the opening of a gap in the particle- 
hole channel at the boundary of the Antiferromagnetic 
Brillouin zone (AFBZ). Note, as argued above, we as- 
sume that this gap can occur in a spin liquid without 
symmetry breaking. We neglect bilayer splitting of the 
bands and work with a single-planar band with disper- 
sion: 

ek = cos kx cos ky — ja 



± Y4t2(cosA:x + cos%) +A|^(k) (1) 

The sign follows the opposite of {coskx + cos ky). 
The gap is assumed to have the form AslO^) = 
Asp [coskx — cos ky -\- a [cos'^ kx -\- cos^ ky]) . If a is 
zero the gap has d-wave symmetry and the same magni- 
tude at (0, tt) and (tt, 0).When a is negative (positive) the 
magnitude of the gap is anisotropic and larger at (7r,0), 
(respectively (0,7r)). The 2a periodicity due to the oxy- 
gen ordering breaks the tetragonal symmetry and splits 
this band into two: a and folding the Brillouin zone 
in the (0, 0) — (tt, 0) direction, and opening a gap V at 
kx = 7r/2, see Bascones et aP^.ln the reduced Brillouin 
zone, kx G (— 7r/2,7r/2) and ky G (— 7r,7r), 
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In the superconducting 



state, the energy of the quasiparticles is E'-^ 
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and the breaking of the tetragonal symmetry the super- 
conducting gap is modified from a pure d-wave form to 



A^'^ = ^ [g{k)coskx t cosky] with 



^(k) = (ek+TT, - Ck) / (ck+TT, - Ck) + V 
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The resonantly enhanced^^ spin susceptibility at low en- 
ergies is given by: 



/mx(q,cj) = JRex'^'''{ci,uj)Imxo{ci,u) 



(4) 



J is a constant of order the exchange energy, /mxo(q, ^) 
is the imaginary part of the bare spin susceptibility of 
the truncated Fermi surface and i?ex^^^(q, cj) is the real 
part arising from the resonant mode. For simplicity, we 
assume that the high energy part of the spin fluctuation 
spectrum is a single mode with dispersion taken from 
experiment^ 



^res(q) = (A 

Here A^^^^=33 meV and A= 
Kronig the real part 
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=365 meV A. From Kramers- 



Rex'-'%q,u;) = l/[Lo- 
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FIG. 1: a) Truncated Fermi surface, plotted in an extended 
Brillouin zone, for the band parameters discussed in the text 
with the tetragonal symmetry only broken by the inclusion of 
the ortho-II ordering, b) Spin susceptibility -convoluted with 
the resolution factor -versus momentum, measured from (tt, tt) 
along ((5, 0) and (0, 8) in solid and dashed lines, respectively, 
for the same parameters and uj — 0.05. The inset shows the 
same plot for the case ty — 1.3tx, see the text 



At zero temperature, applying a Kubo formula. 
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Plus (minus) sign applies when m and m' are equal 
(different), and 



Due to the ortho-II ordering ^k,q 
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Neutron scattering experiments are affected by the reso- 
lution factor of the measurement apparatus. In the ex- 
periments by Stock et a? the resolution factor in q space 
is anisotropic and comparable to the observed incom- 
mensurability. In order to compare with their results 
we convolute ImxiQ^^) with an anisotropic gaussian^^« 
with standard deviation along directions perpendicular 
and parallel to the chains ax = 0.078 and ay = 0.058. We 
set f = 0.15, Asp = 0.31, V = 0.6, and /i = -0.70, in 
units of t, corresponding to hole doping x = 0.09. In the 
following, we assume t = 180 meV, then Aspin = 0.183 
and A = 0.52. 

We first consider the case Aq = a = and analyze 
the possibility that the observed anisotropy is a conse- 
quence of the modification of the Fermi surface due to 
ortho-II ordering. The truncated Fermi surface obtained 
with these parameters is shown, in the extended Brillouin 
zone, in Fig. la). The breaking of the tetragonal symme- 
try is clear. The magnetic response (not shown) has two 
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FIG. 2: a) Fermi surface corresponding to the case which 
includes both saddle-point gap anisotropy and ortho-II order- 
ing (red), and only ortho-II ordering (black). The later one 
is also plotted in Fig. 1. b) Spin susceptibility at a; = 0.05- 
convoluted with the resolution factor- as a function of mo- 
mentum, measured from (7r,7r), along the (^,0) and (0,^), in 
solid and dashed lines, respectively, for a = —0.33 in the nor- 
mal (top curves) and superconducting state with Ao = 0.1 
(bottom curves), c) Colormap of the spin response for the 
same parameters as in b) in the normal state, but without 
including resolution factor effects. 



peaks close to (tt, tt), along the directions (tt + S^tt) and 
(7r,7r + 6). The position of these peaks depends on the 
band parameters and their distance from (tt, tt) increases 
with doping. The spin susceptibility, convoluted with the 
resolution factor, along (tt + J, tt) and (tt, tt -\- S) is shown 
in Fig. lb). As the value of J is uncertain, the suscepti- 
bility is given in arbitrary units, but these units are the 
same throughout. In spite of the clear anisotropy of the 
Fermi surface, once the resolution factor is included, the 
spin response in the directions perpendicular and paral- 
lel to the chains is very similar. We have repeated the 
calculations for different Fermi surface parameters. The 
anisotropy we obtained is very weak or it gives larger 
incommensurability in (tt, tt + (5), contrary to the experi- 
ment. We have checked the effect which, combined with 
the ortho-II ordering, could have the addition of a third 



FIG. 3: Dependence of the spin response -convoluted with 
the resolution factor- on energy for a — —0.33. Top figure: 
color map as a function of energy and momentum, measured 
from (7r,7r), along (^,0) (right figure) and (0,^) (left figure). 
Bottom figure: at (7r,7r). Solid and dashed lines correspond 
to the normal and superconducting states, respectively. 



nearest neighbor hopping term, the reduction of the spec- 
tral weight along the Fermi surface as the saddle points 
are approached or the inclusion of finite lifetime of the 
quasiparticles^"^. None of these effects gives an incom- 
mensurability comparable to the observed value. We con- 
clude that the modification of the Fermi surface due to 
the ortho-II symmetry breaking alone cannot explain the 
experimental results. 

The anisotropy in the hopping matrix element has been 
suggested^^, alone or with other effects, as the origin of 
the one-dimensionality of the spin spectrum reported by 
Mook et a?^ in YBa2Cu306.6- The case in which the 
hopping along {ty) and perpendicular {tx) to the chains 
are inequivalent is shown in the inset of Fig. 1 b). We 
have taken ty = IMx = 1.3t, keeping the rest of the 
parameters as above. With this modification of the Fermi 
surface the tail of the peak at (tt, tt) decays more slowly 
along (7r+(5, tt), than along (tt, 7r-\-S). But this broadening 
is different to the one found in ortho-II YBa2Cu306.5. 
We have also checked the effect of a modulated doping, 
which could follow the ortho-II oxygen ordering, with no 
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success. 

Next we examine the effect of finite a. The gap at 
the AFBZ loses its d-wave symmetry and for a = —0.33, 
at (tt, 0) it is twice than at (0,7r). In spite of this large 
asymmetry, the effect on the shape of the Fermi arc is 
not unreasonable, see Fig. 2 a). The crossing between 
the arcs and the AFBZ moves towards (0, tt) in the first 
and second quadrants and towards (0, — tt) in the third 
and fourth quadrants. The color map in Fig. 2 c) shows 
the spin susceptibility at a; = 0.05, without the reso- 
lution factor. The inequivalence of the peaks is clear. 
This anisotropy remains when the resolution factor is in- 
cluded. The shape of the peaks along (tt + tt) and 
(tt, 7r + (5) are plotted in Fig. 2 b), and compare well with 
those reported^. As shown in Fig. 3 (top figure) the 
anisotropy is maintained with increasing energy. As the 
resonance is approached the incommensurability vanishes 
due to the enhancement factor. The anisotropy is robust 
and survives also in the superconducting state, see Fig. 
2b). In this plot it is clear that at low energies, the mag- 
netic intensity has been reduced in the superconducting 
state. The suppression in the superconducting state is 
due to the opening of a gap along the arcs. It is also 
evident in Fig. 3 (bottom figure) which plots the sus- 
ceptibility at (tt, tt) as a function of energy. The signal 
increases with energy, as reported experimentally. 

In conclusion, we have presented a simple phenomeno- 
logical model which is able to account for the main fea- 
tures of the spin scattering in ortho-II YBa2Cu306.5, 
at low energies, without involving one-dimensionality. 
The model should be applicable to other underdoped 
cuprates. The spin response is larger close to (7r,7r), but 
sensitive to the parameters which determine the Fermi 



surface. It is characterized by two peaks along the Cu- 
O axis. Below Tc the opening of the superconducting 
gap on the arcs produces a reorganization of the spectral 
weight at low energies. The superconducting state can 
influence the spin susceptibility also through a better def- 
inition of the quasiparticles. The model is intended for 
low energies in the underdoped regime. The high energy 
part observed experimentally, including the resonance, 
is modelled with a collective triplet mode . Recently, 
Tranquada et a?^ has suggested that the high energy 
mode is due to disordered bond-stripes, which behave as 
weakly coupled two- leg ladders. However in our view the 
ring form of the dispersion at high energies suggests a 2- 
dimensional character for the spin liquid. Our model can 
be seen as a phenomelogical extension for the pseudogap 
phase of the RVB (resonant valence bond) theory, whose 
success was recently reviewed^^. 

The anisotropy of the spin response found in Ortho-II 
YBa2Cu306.5 is ascribed to an asymmetric deformation 
of the Fermi arcs close to the boundary of the AFBZ. 
We found that such a deformation could originate from 
an anisotropic /S^sl -the energy gap that truncates the 
Fermi surface in the normal pseudogap phase. A sub- 
stantial anisotropy of A^^ was required to explain the 
experiment. The origin of this anisotropy remains to be 
an open question, but could be related to the chain or- 
dering. 
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